function y = spectralDecomposition(sigma)
oxx = sigma(1);
oyy = sigma(2);
oxy = sigma(3);
t1 = (oxx + oyy) / 2;
t2 = sqrt(((oxx - oyy)/2)^2 + oxy^2);
o1 = t1 + t2;
o2 = t1 - t2;
y = [o1 o2];
end
